Phase separation in Bose-Fermi-Fermi Mixtures as a probe of Fermi superfluidity 



o 
o 









I 

o 



en 
> 

o 
i> 

(N 

^' 

O 
OO 

O 



X 



S. G. Bhongale and H. Pu 
Department of Physics & Astronomy, and Rice Quantum Institute, Rice University, Houston, TX 77005, USjQ 

(Dated: June 19, 2008) 

We study the phase diagram of a mixture of Bose-Einstein condensate and a two-component Fermi 
gas. In particular, we identify the regime where the homogeneous system becomes unstable against 
phase separation. We show that, under proper conditions, the phase separation phenomenon can 
be exploited as a robust probe of Fermi superfluid. 
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Mixtures of superfluids open up possibilities of study- 
ing interacting macroscopic quantum systems. The at- 
tempt of such studies started decades ago in the system 
of ^He-^He mixtures. However, the transition into the su- 
perfluid phase of ^He atom in these mixtures occurs at an 
extremely low temperature and has never been reached 
in experiment. Realization of superfluids in atomic quan- 
tum gases makes such studies possible for the first time. 

The atomic analogy of a superfluid '^He-^He system 
is a mixture of Bose-Einstein condensate (EEC) and su- 
perfluid Fermi gas (i.e., a two-component Fermi gas with 
attractive interaction). In this paper, we explore the rich 
phase diagram of this system at zero temperature. By 
studying the free-energy of the homogeneous mixture, we 
identify the regime where the homogeneous mixture be- 
comes unstable against phase separation. Phase separa- 
tion is quite a generic phenomenon occurring in trapped 
atomic mixtures, originating from the interplay between 
the interactions and the spatial variation induced by the 
trap, thus allowing different regions of the trap to fa- 
vor fundamentally different phases [ll, y]. In addition, 
we demonstrate a novel application of the phase sepa- 
ration phenomenon- the detection of Bardeen-Cooper- 
Schricffcr (BCS) superfluidity within the Fermi gas [3|,|j]. 
Such a detection proposal is motivated by the need for 
an efficient superfluidity detection probe for understand- 
ing pairing within unbalanced Fermi mixtures, a scenario 
that has recently attracted a lot of attention. 

The probing concept can be simply understood if we 
envision a BEG localized to a small region within the 
Fermi medium. We assume that the Fermi gas is uncon- 
fined which is a valid approximation if the spatial exten- 
sion of the Fermi cloud is much broader than that of the 
BEG. In the phase separation regime, the BEG may exist 
in a pure form as an isolated bubble surrounded by the 
Bose-Fermi mixture. We show that the formation of such 
a BEG bubble may be made sensitive to the superfluid 
property of the fermionic medium. Thus, the onset of 
BGS superfluidity is signalled by the formation of small 
BEG bubbles which can be readily detected via a den- 
sity measurement. Essentially, bosons serve as a matter- 
wave probe of Fermi superfluid. This idea is illustrated 
schematically in Fig. [TJ 

To begin with, we first determine the zero-temperature 
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FIG. 1: (Color online) Schematic of the proposed BEC probe. 
Bosons are confined in the tight trap and are made to interact 
with the Fermi atoms which are unconfined. Phase separation 
can be easily detected by measuring the BEC density profile. 



projected phase diagram of a Bose-Fermi-Fermi mixture 
comprising of bosons of one species and equal population 
spin-up and -down fcrmions of another. Here, by pro- 
jected we mean a 2-dimensional slice of the d-dimensional 
parameter space by fixing the d — 2 independent param- 
eters. In the absence of inter-fermion interaction, this 
system can be identically mapped onto a Bose-Fermi mix- 
ture and has been previously studied quite extensively [6|] . 
However, the interacting case remains much less explored 
[7| and is of focal importance to this work. Here, we 
therefore assume attractive s-wave interaction between 
fermions of unlike spins. Furthermore, we assume Bose- 
Bose and Bose-Fermi interactions to be repulsive. It is 
convenient to treat all interactions via a pseudo-potential 
modeled by Va{r — x) = AQ,i5(r — x) which is valid for 
dilute systems, where we indicate different types of in- 
teracting atoms by the index a € {BB,B 1,B i,ti}- 
In this work, we make two non-essential simplifications. 
First, we assume that the Bose-Fermi interaction is spin- 
independent, i.e., Ab| = Ab| = Xbf- Second, we con- 
sider a single spatial dimension. We remark that our 
theoretical framework is general and the above two re- 
strictions can be straightforwardly removed. 

The homogeneous system under study is characterized 
by a total of 5 independent parameters: the 3 interac- 
tion strengths {Xbb, Xbf, X-h) and two densities (ps, gp) 
where g^ — Qi — Qf/2. Note that the choice is not 
unique. For instance, instead of the densities, we may 
choose chemical potentials {p,b,^.f) to be independent 



parameters. It is instructive to first consider the much 
simpler case without Fermi-Fermi interaction, or A|| — 0. 
As mentioned earher, this case, for 3D confinement, has 
been studied by several authors [6|. Similar analysis can 
be performed for the ID case, resulting in the instabil- 
ity criterion given by Qp < '^\%p/{tt'^Xbb) — when this 
inequality is satisfied, the homogeneous mixture is unsta- 
ble and tends to phase separate. Thus the BEC density 
profile in the two regimes (phase separated or not) may 
be made quite distinct as depicted in the schematic of 
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Fig. [T] However, we emphasize that, since the interac- 
tion is density-density, there is no direct connection of the 
fermionic pairing gap (in the case of interacting fermions) 
to the phase- separation phenomenon. Remarkably, as we 
shall show, the phase separation can be made sensitive 
to the fermionic pairing, underlying the basis of our pro- 
posed BEC probe of Fermi superfluidity. 

We begin by writing the free energy functional of the 
homogeneous mixture in the form, 
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where we have defined ek = h, fc^/2m, jlp = t^F — ^BFQB, 
and the fermionic quasi-particle energies are given by 
Afe = ^(£fe — JIfY + ^■^ by neglecting the Hartree con- 
tribution since it only leads to a constant energy shift. 
Also, anticipating the role of Fermi supcrfluid, we have 
written the free energy as a function of boson density qb 
and BCS pairing gap A. We will then construct the pro- 
jected phase diagram in the qb-^ parameter space while 
fixing the values of three more independent variables to 
be discussed below. 

The thermodynamic ground state is given by the min- 
imum of the the free energy and therefore corresponds 
to the necessary first derivative conditions: 9F/9A = 
and dF/dgB = 0. The first of these conditions 
essentially reproduces the gap equation in the form 
(— A|j^/27r) /g°°(l/Afe)dA; — 1, while the second fixes the 
number through the modified bosonic Thomas-Fermi 
equation XbQb — fJ'B + XbfQf = 0. However, local min- 
imum is guaranteed only if the Hessian matrix M, con- 
structed from the second derivatives of F, is positive def- 
inite, or the following conditions are satisfied: 



M„ 



a^F 

Wb 



>0, Det[M]>0. 



(2) 



When condition ^ is violated, the system will neces- 
sarily phase separate [8|. Using this simple criteria, one 
can map the whole phase space of the homogeneous mix- 
ture. However, this is an extremely laborious task, given 
that the total phase space is huge, represented by the set 
of five independent parameters as we mentioned earlier. 
We make a judicious choice and use the following set: 
{Abb, Abf, A, Qb, ^J-f}, the motivation behind which will 
be clear as we proceed further. Other parameters, such 
as A-fi, gp and hb, must be calculated self-consistently 
using the gap equation and the bosonic Thomas-Fermi 
equation discussed earlier, together with the fermionic 
number equation gp = (2/7r) J„ [1 — {sk — jlp)/ IS-t^]dk. 



We remind the reader that our goal here is not to give 
a complete description of the whole phase space but focus 
our study on a small region that is physically meaning- 
ful in view of current experimental setups and illustrate 
a novel probing techniques based on the phase separa- 
tion phenomenon. This we do by first picking reasonable 
values of Xbb^ Xbf, and iip. We then determine the 
stable/unstable regions in the ps-A space via condition 
(121) . For this we need to study the properties of the Hes- 
sian matrix M. States that satisfy condition ([2]) are only 
guaranteed to be a local minimum of the free energy. To 
determine whether the state is the ground state of the 
system, we need to compare the free energies of different 
homogeneous phases which include the pure BEC and 
the pure Fermi phase, in addition to their mixture. 

Following the above procedure, we have determined 
the phase space for various values of the fixed parame- 
ters. We find that this system exhibits a very rich phase 
diagram. However given the lack of space, we restrict 
ourselves to pointing out some general features that are 
relevant to this proposal and direct the reader to an up- 
coming publication for more details. 

The projected phase diagram of a homogeneous mix- 
ture for a particular set of {Xbf,Xbb, ^J'F} is shown in 
Fig.[2fa). We note the following: (1) The unstable region 
correspond to a partial ellipse cut by the gB axis (shown 
by the shaded region in Fig.[2Ka)). (2) Within this region, 
the phase space can be further divided into a dynamically 
unstable region (the yellow/brighter shaded) at which the 
free energy landscape shows a saddle point and a dynam- 
ically stable region (the green/darker shaded) at which 
the free energy landscape shows a local but not a global 
minimum. (3) The position and the extent of these re- 
gions depends on the values of the fixed parameters. Im- 
mediately we notice that, in a typical experimental setup 
with a given value of A||, only a small part of the phase 
space comprising of points on the fixed A|j contour, for 
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FIG. 2: (Color online) (a) Projected phase space of the Bose-Fermi-Fermi mixture. Here, as well as in other figures, all 
quantities are scaled in units of the probe, represented by an harmonic oscillator with frequency loq and length Iq. Fixed 
parameters are \bf = 0.44, \bb = 0.05, and fip ~ 27.32. The yellow (green) area indicates region where the Bose-Fermi 
mixture is dynamically (mechanically) unstable. The contour, Ci and C2, corresponding to A|| = —4 and —6 respectively, 
(b) Comparison of free energy of the mixture (solid curve) along the contour Ci and the free energy of the pure BEC (dashed 
curve). The free energy the pure Fermi gas is much higher and not shown. The points shown by red symbols correspond to 
those in (a), (c) Critical boson chemical potential, p,B as a function of the inter- Fermi interaction separating the pure BEC and 
mixed Bose-Fermi phase. Within LDA, moving vertically upwards along a line parallel to the y axis implies moving from the 
edge towards the center of the trap. 



example Ci or C2 shown in Fig.[5Ja) is physically acces- 
sible. The vertical axis of A = corresponds to a mixture 
of BEC with a non-interacting Fermi gas (Af| = 0). The 
significance of the particular parameter set chosen ear- 
lier is clear since we can now directly obtain the Bose 
density profile in the probe by mapping the boson chem- 
ical potential along this contour, /^^[Ci] on to the spatial 
coordinate in the probe via the local density approxima- 
tion (LDA) using fiB[Ci] — fJ-sir) = //b — V"(r[Ci]) where 
V{r) is the probe trapping potential for BEC. In Fig.[2I^b) 
we plot the free energy of the mixture as a function of 
the Bose chemical potential ij,b[Ci]. The same plot also 
show the free energy for a pure BEC. Now we can eas- 
ily identify the following special points in Fig. [2l[a): the 
point marked by the circle (triangle) as the point where 
the free energy of the BEC goes below that of the mix- 
ture and thus represents a first order phase transition in 
the presence(absence) of pairing, and the point marked 
by square as the point where hb reaches its maximum 
value along the contour Ci. 

Now we are in the position to discuss how we can take 
advantage of the phase diagram to detect Fermi super- 
fluidity using BEC as a probe. To put this idea in con- 
text, we note that there is a new avenue in cold atom 
research that cast BEC as tools for quantum measure- 
ment. A good example is the experiment by Kriiger et 
al. [9[ where the Thomas-Fermi character of the BEC 
density profile is exploited to measure the surface po- 



tential energy landscape with exquisite accuracy. Also, 
very recently, Bhongale and Timmermans have proposed 
a high sensitivity force measurement by exploiting phase 
separated two-BEC mixture [10| . The ideas in this letter 
are a next step in this direction. 

Since for our system the bosons are confined to the 
probe, we want a situation such that, on phase separa- 
tion, the system consists of a pure BEC component near 
the center of the trap surrounded by a cloud of either 
mixed bosons and fermions or pure fermions. It is quite 
intuitive that if //^(r) > /iB[A|J (the latter being the 
boson chemical potential along the stable boundary of 
the green region represented by the circle in Figs. [H^a) 
and [DJb) , a pure BEC bubble will phase separate out 
from the mixture. However, in order to use this phe- 
nomenon to discern the fermion superfluid property, it is 
important to connect the occurrence of the BEC bubble 
with the disappearance/ appearance of a nonzero super- 
fluid gap A. This is crucial since, as mentioned earlier, 
phase separation can also occur in a mixture of BEC and 
a normal Fermi gas due to Bose-Fermi repulsion. For this 
reason, in Fig. ^c) , we plot the critical boson chemical 
potential as a function of the inter-fermion interaction 
strength. We see that there is a clear separation of re- 
gions corresponding to a pure BEC phase above the curve 
and a mixed phase below. Thus if we configure the probe 
such that the center chemical potential is very close to 
and above the curve in Fig. ^c) implying a phase sepa- 




FIG. 3: (Color online) Boson density profile within LDA with 
(dashed) and without (solid) fermion pairing, for total number 
of bosons Nb = 101, at A|| = —A. All other parameters are 
same as used in the previous plot. 



rated pure BEC component near the center of the trap, 
the appearance of BCS superfluidity on increasing the 
attractive inter-ferniion interaction is immediately sig- 
naUed by the disappearance of the phase separated pure 
BEC bubble as shown in Fig. [3] 

In, practice, the sharp density jump will be smoothed 
out due to the finite kinetic energy contribution, however, 
the absence of non-Thomas Fermi variation in the Boson 
density profile can be easily detected. Also, we claim that 
this probing method allows for a numeric estimation of 
the gap. A, by a systematic fitting technique if the exact 
density profile is obtained via a numerical solution of the 
coupled Bose-Fermi equations. 

Moreover, increasing the attractive interaction beyond 
a certain value (depending on the value of Xbb and \bf) 
results in the absence of a phase separated regime as de- 
picted by the contour C2 in Fig. [Ifa). In fact, this can 
be considered as a very strong indication of the pairing 
phenomenon, the limiting case of which is a homogeneous 
mixture of molecular (M) and atomic BEC on the repul- 
sive side of the Feshbach resonance, known to be com- 
pletely stable if \bm < V^bb^mm [ll- However, this 
involves calculating these additional interaction strength, 
which will be dealt with in a future article. 

In conclusion we have studied the phase diagram of 
a Bose-Fermi-Fermi mixture and proposed a probe for 
detecting BCS type superfluidity within a quasi-lD two- 
component Fermi gas by configuring the system such that 
thermodynamic instability resulting in phase-separation 
is sensitive to the interaction between the fermions and 
hence the BCS pairing. The probe consists of a BEC con- 
fined to a relatively tight trap. We have shown that by 
properly tuning probe parameters, the density profile of 
the BEC provides a robust signal of the fermion pairing. 
The probe idea may be easily extended to 3D systems 
by identifying appropriate phase separation regime. The 
expressions for the free energy in Eq. ([1]) and the cor- 
responding Hessian matrix are still valid in 3D with the 
understanding that the integrals involved are also 3D. 
One important difference between 3D and ID is that, in 
the former, proper renormalization procedures must be 



taken to remove the ultraviolet divergence in gap equa- 
tion associated with the contact interaction. 

We are aware of a related recent proposal for probing 
BCS type superfluidity using a overlapping BEC. How- 
ever there, the pairing signal is related to the damping 
of the BEC acoustic phonons ll|. The strength of our 



proposal lies in the fact that the signature of pairing is 
reflected in the BEC density profile, a quantity that can 
be easily measured in experiment. The proposed prob- 
ing scheme possesses another important advantage: it 
probes the local value of A. This is crucial in trapped 
experiment where the gap varies in space due to the trap- 
induced inhomogeneity, and in situation where the gap 
has intrinsically nontrivial spatial dependence. The lat- 
ter arises, for instance, in the case of FFLO superfluid 
state in a population imbalanced Fermi system where 
the gap varies sinusoidally in space, which may result 
in density oscillations in BEC probe. Another natural 
extension is to study the phase diagram of the Bose- 
Fermi-Fermi mixture where the Fermi-Fermi interaction 
is tuned across a Feshbach resonance. How the presence 
of the bosons affect and probe the BEC-BCS crossover 
will be an interesting problem to study. 
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